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The relative isobaric yields of fragments produced in a series of heavy ion induced multi- 
fragmentation reactions have been analyzed in the framework of a Modified Fisher Model, primarily 
to determine the ratio of the symmetry energy coefficient to the temperature, Oa/T, as a function of 
fragment mass A. The extracted values increase from 5 to '-^ 16 as A increases from 9 to 37. These 
values have been compared to the results of calculations using the Antisymmetrized Molecular Dy- 
namics (AMD) model together with the statistical decay code Gemini. The calculated ratios are 
in good agreement with those extracted from the experiment. In contrast, the ratios determined 
from fitting the primary fragment distributions from the AMD model calculation are ~ 4 and show 
little variation with A. This observation indicates that the value of the symmetry energy coefficient 
derived from final fragment observables may be significantly different than the actual value at the 
tim e of fragment formation. The experimentally observed pairing effect is also studied within the 
same simulations. The Coulomb coefficient is also discussed. 
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I. INTRODUCTION 

In the early 1980's, a study of the isotopic yield distri- 
butions of intermediate mass fragments produced in high 
energy proton induced multi-fragmentation reactions at 
Fermi Lab showed that the distributions can be weh de- 
scribed by a Modified Fisher Model [li,|2|] in which the iso- 
tope production is governed by the available free energy. 
Therefore isotope yields provide a good probe to study 
the nature of the disassembling nuclear system. Multi- 
fragmentation of the system is also generally observed 
in violent collisions in heavy ion reactions in the Fermi 
energy domain and there is evidence that both subnor- 
mal and supernormal densities may be explored in such 
collisions 0, 0] • Work in this area has concentrated on 
exploring the nuclear equation of state and the liquid gas 
phase transition in nuclear matter. In a recent paper we 
addressed the possibility of probing the quantum nature 
of the liquid-gas phase transition using isotope distributi 
ons 0. 

Over the last several years many fragment emission 
studies have been motivated by efforts to use fragment 
yield distributions, either singly or by comparison to 
those of similar reactions, to explore the symmetry en- 
ergy in the emitting source at different densities and tem- 
peratures d, . In each of these cases measuring the 
isotope distributions over a wide range of mass number A 
and atomic number Z should provide a more reliable basis 



for extraction of the desired information. Even then there 
are important issues which must be resolved in order to 
establish the relation between the experimental isotope 
distributions and the symmetry energy of the emitting 
system. One is the source temperature T. In the Modified 
Fisher Model as well as other approaches based upon the 
free energy, all terms which can be determined from ex- 
periments appear in the form of /T times some function 
of A, Z or the neutron number N, where i indicates the 
coefficient of the different terms contributing to the free 
energy. Since the beginning of the experimental study 
of heavy ion collisions in the multifragmentation regime, 
significant efforts have been made to evaluate the source 
temperature, but no absolute consensus among different 
methods has yet been achieved A second issue is 

the role and effect of secondary decay processes. In ex- 
periments the majority of the detected fragments are in 
their ground states. Most of the primary fragments pro- 
duced in Fermi energy heavy ion reactions are expected 
to be in an excited state when they are formed. Indeed, 
in previous work, excitation energies of the primary frag- 
ments have been evaluated by studying the associated 
light charged particle multiplicities ll|, |l2[. Such data 
demonstrate that secondary decay is important and raise 
the question of the degree of confidence which can be ac- 
corded to experimental derivations of the symmetry en- 
ergy coef ficient which do not properly correct for this 
important effect. A clear goal for experimentalists would 
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be the reconstruction of the primary isotope distribu- 
tions from the experimental distributions. This might 
be approached by a reconstruction of the primary iso- 
tope distributions employing the associated neutron and 
charged particle multiplicities. However, since multiple 
fragments are produced in a reaction and light particles 
can be produced even before the formation of the frag- 
ments, the identification of the parent for a detected light 
particle observed in coincidence with detected fragments 
is not straight forward. A reconstruction of the primary 
isotope distribution was part of goal of the experiment 
described here, but that analysis is still underway [l^l ■ 

A third issue, not addressed in this paper but of critical 
importance in this field, is the problem of obtaining reli- 
able estimates of the density at the time of fragment for- 
mation. In the absence of such information most density 
estimates should be viewed as unconfirmed. While this 
may be somewhat mitigated by comparing experimental 
observables with results of dynamic models employing 
particular assumed forms for the density dependence of 
the symmetry energy, such approaches are integral and 
may be influenced by other assumptions and parameter 
choices inherent in the model applied 14 1. 

In this paper we explore the extent to which informa- 
tion on the symmetry energy, in the form of aa/T, can 
be extracted from high quality data for isotope resolved 
fragment yield distributions and compared to the model 
predictions The role of the secondary decay is 

explored by comparisons with results of theoretical cal- 
culations. In a future paper we will discuss the extraction 
of such information using isoscaling techniques 16[. 



II. EXPERIMENT 

The experiment was performed at the K-500 super- 
conducting cyclotron facility at Texas A&M University. 
^''^™Zn and ^''Ni beams were used to irradiate ^^'^^Ni, 
"2,i24sn, and 232Th targets at 40 A MeV. Inter- 

mediate mass fragments (IMFs) were detected by a de- 
tector telescope placed at 20°. The telescope consisted 
of four Si detectors. Each Si detector was 5cm x 5cm. 
The nominal thicknesses were 129, 300, 1000, 1000 /im. 
All Si detectors were segmented into four sections and 
each quadrant had a 5° opening angle in polar and az- 
imuthal angles. Therefore the energies of the fragments 
were measured at two polar angles of the quadrant de- 
tector, namely 9 = 17.5° ± 2.5° and 9 = 22.5° ± 2.5°. 
Typically 6-8 isotopes for atomic numbers, Z, up to Z=18 
were clearly identified with the energy threshold of 4-10 A 
MeV, using the AE-E technique for any two consecutive 
detectors. The De^iaE-E spectrum was linearized em- 
pirically. Mass identification of the isotopes were made 
using a range-energy table [l3| . In the analysis code, iso- 
topes are identified by a parameter Z^eai ■ For the isotope 
with A=2Z, Zneai = Z is assigned and other isotopes are 



identified by interpolation between them. Typical Zneai 
spectra are shown in FiglT] The energy spectrum of each 
isotope was extracted by gating the isotope in a 2D plot 
of Zneai vs energy. The yields of light charged particles 
(LCPs) in coincidence with IMFs were also measured us- 
ing 16 single crystal CsI(Tl) detectors of 3cm thickness 
set around the target. The light output from each de- 
tector was read by a photo multiplier tube. The pulse 
shape discrimination method was used to identify p, d, 
t, h and a particles. The energy calibration for these 
particles were performed using Si detectors (50 -300 /im) 
in front of the Csl detectors in separate runs. 

The yield of each isotope was evaluated, using a moving 
source fit. For LCPs, three sources (projectile-like(PLF), 
nucleon-nucleon-like(NN) and target-like (TLF)) were 
used. The NN-like sources have source velocities of about 
a half of the beam velocity. The parameters are searched 
globally for all 16 angles. For IMFs, since the energy 
spectra were measured only at the two angles of the quad- 
rant detector, the spectra were parameterized using a 
single NN-source. Using a source with a smeared source 
velocity around half the beam velocity, the fitting param- 
eters were first determined from the spectrum summed 
over all isotopes for a given Z, assuming A=2Z. Then 
all extracted parameters except for the normalizing yield 
parameter were used for the individual isotopes. This 
procedure was based on the assumption that, when the 
spectrum is plotted in energy per nucleon, the shape of 
the energy spectrum is same for all isotopes for a given 
Z. Indeed the observed energy spectra of isot opes are 
well reproduced by this method. For IMFs, a further 
correction was made for the background. For the IMFs, 
a further correction was made for the background. As 
seen in FiglTJ the isotopes away from the stability line, 
such as and ^^P, have a very small yields and the 
background contribution is significant. In order to eval- 
uate the background contribution to the extracted yield 
from the source fit, a two Gaussian fit to each isotope 
combined with a linear background was used. The fits 
are shown in FiglTJ Each peak consists of two Gaussians. 
The second Gaussian (about 10% of the height of the 
first one) is added to reproduce the shape of the valley 
between two isotopes. This component is attributed to 
the reactions of the isotope in the Si detector. The cen- 
troid of the Gaussians was set to the value calculated 
from the range-energy table within a small margin. The 
final yield of an isotope with Z > 2 was determined by 
correcting the yield evaluated from the moving source fit 
by the ratio between the two Gaussian yields and the 
linear background. Rather large errors ( ^ ±10%) are 
assigned for the multiplicity of the NN source for IMFs, 
originating from the source fit besides the background 
estimation. The errors from the source fit are evaluated 
from the different assumptions of the parameter set for 
the source velocity and temperature. 
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FIG. 1: Typical linearized isotope spectra for Z=6 and 15 are 
shown for ^*Ni+^-^*Sn. The number at the top of each peak is 
the assigned mass number. The histograms depict experimental 
data. A linear back ground is assumed from valley to valley for a 
given Z. Each isotope is fit with two Gaussians. The individual fit 
indicates the yield of the isotope above the back ground. The sum 
of Gaussians and the background are also shown in each spectrum. 



III. MODIFIED FISHER MODEL 

In order to study the symmetry energy contribution 
to the isotope production, the Modified Fisher Model of 
ref.l is used. In the Modified Fisher model, the fragment 
yield of A nucleons with I=N-Z, Y(A,I) is given by 



Y{A,I) = CA-^exp{[{W{A,I) + finN - 
+Nln{N/A) + Zln{Z/A)}. 



f,,Z)/T] 



C is a constant. The A""^ term originates from the en- 
tropy of the fragment and the last two terms are from the 
entropy contributions for the mixing of two substances in 
the Fisher Droplet Model [l^]- is the neutron chem- 
ical potential and is the proton chemical potential. 
W(A,I) is the free energy of the cluster at temperature 
T. As such it includes both energy and entropy terms. In 
the model W(A,I) is given by the following gene ralized 
Weiszcicker-Beth semi-classical mass formula 1^, 2^ at 
a given temperature T and density p, 



WiA,I) 



-aa{p, T)f/A ~ a,ip, T)Z{Z - l)IA^/^ 
-a,{p, T)A - as{p, T)A'/^ - 5{N, Z). (2) 



(1) 



The indexes, v, s, c and a represent volume, surface. 
Coulomb, and symmetry energy, respectively. Following 
the semi-empirical mass formulation, the pairing energy, 
(5(N,Z), is given by i21,,i22,] 



ap(p,r)/^i/2 (odd-odd) 
5{N, Z)^ { (even-odd) (3) 

- ap[p,T)/A^/'^ (even-even). 



We define the isotope yield ratio, R(I-|-2,I,A), between 
isobars differing by 2 units in I as 



R{I + 2, /, A) ^ Y{A, I + 2)/Y{A, I) 
= exp{[{W[I + 2, A) - W{I, A) + {iin - tip)yT + 5„,„(/ + 2, A) - 5„„(/, A)}, (4) 



where Smix{I, A) = Nln{N / A)+ Zln{Z / A) . Hereafter, in as ai= a.i{p,T) (i=v,s,c,a,p). Inserting Eq.Q into Eq.©, 
order to simplify the description, the density and temper- one can get 
ature dependence of the coefficients in Eq.© is omitted 



R{I + 2, /, A) = exp{\^in - Mp + 2flc(^ - l)/A^'^ ~ Aaa{I+l)/A -6{N^\,Z-1)- 5{N, Z)]/r+ A(/ + 2, /, A)}, (5) 

I 



where A(/ -t- 2, /, A) = S'„„(/ + 2, A)- 5™„(/, A). One 
should note that A(l, — 1, A) = and for other I values 
A(/ + 2,1, A) < 0.5, which is rather small comparing 



other parameters in Eq.(l5]). 

Initially we focus on the isobars with I=-l and 1. For 
these isobars the contributions from the symmetry term 
and the mixing entropy term in Eq.(IS]) drop out and, 
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FIG. 2: Experimental average values of ln[R(I+2,I,A)] for the case 
of I=-l are plotted as a function of Z/A of the reaction systems. 
Z/A = (Zp+Z()/(Ap+Ai) where p and t represent the projectile 
and the target, respectively. The dotted line is a linear fit. 

since these isobars are even-odd nuclei, the pairing term 
also drops out. Taking the logarithm of the resultant 
equation, one can get 

ln[R{l, -1, A)] = [(^„ - Mp) + 2a,{Z ~ l)/A'/'']/T. (6) 

For different reaction systems, N/Z and thus {fin — fJ-p)/T 
can be different. In order to evaluate the system depen- 
dence of {fin — fJ'p)/T we determined the average value 
of ln[R(l,-l,A)] over all available fragments for each sys- 
tem. In Figure [5] these average values are plotted as a 
function of the entrance channel Z/A for the reaction sys- 
tems studied. As seen in the figure, the average values 
show a linear dependence on the entrance channel Z /A of 
the system. Since the Coulomb energy in the right hand 
side of Eq.® is that of the fragment itself and there- 
fore expected to be very similar for the different reaction 
systems and the temperature T is also assumed similar 
because the same incident energy is used, we attribute 
the linear dependence seen in the figure to the differ- 
ence of {fin — ^J■p)/T in the different systems. Expressing 
(Mn ^ ^J'p)/T as 

{[fin ~ fip)/T], = - fXp)/T]o + Ati{Z/A)/T, 

and 

Afi{Z/A)/T = cl • {Z/A) + c2, (7) 

here [fin — ^J'p)/T]i denotes the value for a given re- 
action system i, and — fj,p)/T]o is for the refer- 
ence reaction system. A linear fit to the expression ([7]) 
gives cl=-13.0 and c2=8.7 for Fig.©. We took the 
^^Zn-|-^^^Sn reaction as the reference, i.e., the extracted 
values in the figure have been adjusted to the reference 
reaction using Afi{Z/A)/T, in which Afi{Z/A)/T^O for 



FIG. 3: Experimental values of ln[R(l,-l,A)] with the offset cor- 
rection for different reactions is plotted as a function of A for I=-l. 
Open circles show results from the individual experiments and solid 
circles depict the average values for a given A over all reactions. 
The dotted line shows the result of fitting the average values with 
Eq.®. 

the ^''Zn-l-^^^Sn reaction. In Figl3]the experimental val- 
ues of ln[R(l,-l,A)] corrected by A^(Z/A)/T are plotted 
for all reactions as a function of A. Fitting these cor- 
rected average values using — /ip)/T and ac/T as fit- 
ting parameters, Eq.® leads to (/x„ — /ip)/T— 0.71 and 
ac/r=0.35. 

We next compare isobars with 1=1 and 3, noting that 
these isobars are also even-odd nuclei for which the pair- 
ing term is 0. For this combination, the symmetry energy 
coefficient term in Eq.® is given as a function of A by 

aa/T = -A/8{ln[R{3, 1,A)] - [{fin - fip)/T 
+2a,{Z - l)/A^/^]/T - Afi{Z/A)/T - A(3, 1, A)}. (8) 

In FigH] values of Ua/T calculated from Eq.® using 
the values (/i„ — /ip)/r=0.71 and ac/T—0.35 determined 
above, are plotted as a function of A. All available values 
from the different reactions are plotted in the figure. The 
extracted values are very similar in magnitude and trend 
for the different reactions. In general the values increase 
from 5 to ~ 16 as A increases from 9 to 27 and may show 
a plateauing above that. 

The symmetry term can also be extracted without 
evaluating the values of — fip)/T and Uc/T explic- 
itly. In FiglSl the experimental values of ln[R(3,l,A)] and 
ln[R(l,-l,A)] from the ^''Zn-K^^gn reaction are plotted. 
The symmetry term, aa/T ^ for a give A can be extracted 
approximately by the difference of these values as 

aa/T A/%{ln[R{-i, 1, A)] - ln[R{l, -1, A)] 

-A(3,1,A)}. (8') 

The approximation made in Eq.(8') is that the Coulomb 
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FIG. 4: Experimental values of aa/T as a function of A. Open 
circles are obtained from Eq. (|8ll and solid circles are the average 
values for a given A. Stars are the average values obtained from 
Eq.(8'). 
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FIG. 6: Extracted values of (sign)ap/T are plotted as a function 
of A. Open circles are obtained, using Eq. (|9ll on the upper and 
Eg. moi l on the lower, respectively, for individual reactions. Solid 
circles are averaged values for a given A over all reactions. Stars 
are obtained from Eqs.JTTJ and II12I I. 



FIG. 5: Experimental values of ln[R(H-2.I,A)] for I=-l(opcn cir- 
cles) and I=l(solid circles) for the ^''Zn-f-'^-'^^Sn reaction. 



term in ln[R(3,l,A)] is same as that in ln[R(l,-l,A)]. In 
the actual calculation, the Coulomb term in ln[R(3,l,A)] 
for A=13, for example, is calculated from the yield ratio 
of ^^B/^^C, whereas that in ln[R(l,-l,A)] is calculated 
from ^^B/^^C, assuming that the ratio of the Coulomb 
energy for ^'^B/^'^C is same as that of ^^B/^^C. A similar 
approximation has been made for all extracted values. 
The resultant symmetry values are plotted in Fig fusing 
star symbols. The main difference between the values 
from Eq.® (circles) and those from Eq.(8') (stars) orig- 
inates from the deviation of the data from the fitted line 
in FigH 

The pairing terms in Eq.® can be determined using 
the extracted values for (/i„ — fip)/T, ac/T and aa/T for 



the combination of isobars with 1=0 and 2 and with 1=2 
and 4. For 1=0 and 2 isobars, the pairing term can be 
written as 

Op/T = {sign){l/2)A^'^{ln[R{2, 0, A)] - [(^„ - 

+2a,(Z - 1)7^1/3 - AaalA]/T - A(2, 0, A)}, (9) 

and for 1=2 and 4 it is given by 

Op/T = {sign){ll2)A^/^{ln[R{A, 2, A)] - [(^„ - ^Xp) 

+2ac{Z - l)/A^/^ - 12aa/A]/T - A(4, 2, A)}. (10) 

Here sign=l for (N,Z) = (odd,odd) and -1 for (even, even) 
nucleus, in which A=N-|-Z. The ap/T values obtained 
from Eqs.Q and ([TOl) . using the extracted values of (^n — 
^p)/T, Uc/T and Oq/T above, are plotted in FiglHl On 
the top for 1=0 and 2 isobars, the pairing contribution 
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is clearly observed for isobars with A < 30, though the 
even-odd oscillation pattern is slightly distorted. On the 
other hand on the bottom, using Eq. p^ with 1=2 and 4 
isobars, only clear pairing effect is observed for isotopes 
only with A < 20. 



J 



The pairing term, ap/T, can also be extracted from 
experimental yield ratios of isobars without the explicit 
evaluation of (/i„ — /ip)/T, a^T and aa/T, similar to 
Eq.(8'). Inserting Eq.® and (8') into Eq.® with 1=0 
and 2 isobars, one can get 



flp/T - {sign){l/2)A^/'^{ln[R{2, 0, A)] - ln[R{l, -1, A)] - l/2(?n[i?(3, 1, A)] - hi[R{l, -1, A)] 

-A(3,1,A))-A(2,0,A)} 
= {sign){l/2)A^/^{ln[R{2, 0, A)] - l/2{ln[R{l, -1, A)] + ln[R{i, 1, A)] 

-A(3,1,A))-A(2,0,A)}. 



From Eg. dTUl) with 1=2 and 4 isobars. 



(11) 



ap/T - {sign){l/2)A^/'^{ln[R{A, 2, A)] - ln[R{l, -1, A)] - 3/2(?n[i?(3, 1, A)] - ln[R{l, -1, A)] 

-A(3,1,A))-A(4,2,A)} 
= {sign){l/2:)A^/^{ln[R{A, 2, A)] + l/2{ln[R{l, -1, A)] - 3/n[i?(3, 1, A)] 

+3A(3,1,A))-A(4,2,A)}. 

I 



(12) 



The resultant values are plotted in FiglS] using star sym- 
bols. These are the averaged values over all reactions. 
The results are consistent to those from Eqs.® and ([T0|. 



IV. COMPARISONS WITH MODEL 
CALCULATIONS 

In the multi-fragmentation regime of heavy ion reac- 
tions, fragments may be formed in excited states lll.ll2|. 
Such fragments will de-excite by statistical decay pro- 
cesses. The experimentally detected fragments are nor- 
mally in the ground state. In order to study the ef- 
fect of the secondary decay process on the experimen- 
tally extracted symmetry energy coefficient, we have used 
an Antisymmetrized Molecular Dynamics (AMD) Code 
0, m, 24 [ to model the reaction dynamics and coupled 
it with the statistical decay code Gemini [2^ to model 
the secondary decay processes. The AMD code has pre- 
viously been used to study the fragment production in 
Fermi energy heavy ion reactions and the global fea- 
tures of the experimental results have been well repro- 
duced 26143 1| . We believe that the dynamics can be cru- 



cial in accounting for fragment production in early stages 
of the reaction and the use of a dynamic model , such as 
AMD is essential. 

Since the AMD calculation requires a lot of CPU 



time, only two of the experimental reaction systems have 
been studied in detail. The systems examined were 
64Zn+"2Sn and 64Ni+i24sjj^ ^q^j^ 49 a MeV. Ah re- 
sults shown in this paper have been calculated using a 
newly installed computer cluster in the Cyclotron Insti- 
tute [32[. The calculations have been performed, using 
the Gogny interaction with an asymptotic stiff symme- 
try energy term [9|] , although very similar results are ob- 
tained for the standard Gogny interaction. In order to 
obtain yields of the final products, the yields of primary 
fragments are first evaluated at a given time in the offline 
analysis. The fragments are formed using a coalescence 
radius, i?c, in phase space. Because the hot and dense 
composite system formed at an early stage of the reaction 
expands very quickly, the primary fragment distributions 
are rather sensitive to the choice of coalescence radius and 
the tim e of its application When a smaller R^ is used, 
one can form fragments at an earlier stage. In previous 
calculations, Rc= 5 and t=300fm/c were used and the 



27, m 



experimental results were well reproduced 
Rc—5 corresponds to a radius of 5 fm in configuration 
space. In order to study the effect of the choice of these 
parameters, two different coalescence radii, Rc = 1.5 and 
5, are used here. In the case of Rc=5, the fragment for- 
mation is evaluated at t = 300fm/c. For Rc = 1.5, the 
evaluation is at t = 150fm/c. The excitation energy of a 
fragment is calculated by subtracting the binding energy 
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FIG. 7: (Upper) Calculated multiplicity distributions of the pri- 
mary fragments evaluated using _Rc=l-5 and 5. See details in the 
text. (Lower) Multiplicity distributions of fragments for the exper- 
iments and the calculations. The experimental values are shown 
by open circles. The AMD+Gemini calculated values filtered by 
the experimental acceptance are shown for iJc=1.5 (squares) and 
i?c=5 (triangles). All errors evaluated are smaller than the size of 
the symbols. 



from the total energy. For each isotope the binding en- 
ergy is calculated within the AMD code using a stochas- 
tic cooling method js^l . In the upper part of Figlll the 
primary fragment distributions are shown as a function 
of the fragment Z for the two different cases. As one 
can see, for Rc = 1.5, the multip licity of light IMFs 
with Z < 10 is significantly enhanced, compared to that 
for i?c= 5.0, whereas the heavier fragment yields are sup- 
pressed as expected. The deexcitation of these primary 



fragments was followed using the Gemini code [25|. The 
Gemini code has been used extensively with the AMD 
simulations in the past and goo d ag reement with the ex- 
perimental data has been seen 22, IS 3] • In order to 
sample all possible decay channels of the excited frag- 
ments, one AMD event is used 100 times in Gemini, with 
different random number seeds. In the lower part of Fig[7] 
the multiplicity distributions of the secondary fragments 
are compared with the experimental values. Since the 
experimental values are taken from the NN-source com- 
ponent of the energy spectra, the calculated energy spec- 
tra are subjected to the same filter. For simplicity the 
calculated range of the impact parameter is set to < 
b < 8 fm to suppress the contribut ion from the pro- 
jectile like fragments from the peripheral collisions. The 
impact parameter range was determined from the cor- 
relation between the collision centrality and impact pa- 
rameter studied in ref.30. The angle range, determined 
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FIG. 8: ln[R(I+2,I,A)] for I=-l and 1 primary (open circles) and 
secondary (filled squares) fragments of AMD-Gemini events for the 
^"'Zn-l-li^gjj reaction. A;i(Z/A)/T=0 for this reaction. Solid cir- 
cles are the experimental values for this reaction. Dotted, dashed 
and Dashed-dotted lines show the fits, using Eq. JB]! for the experi- 
ments, the primary and the secondary fragments, respectively. 



from the extracted moving source parameters, was set to 
b° < 9 < 45° to suppress the heavy projectile-like con- 
tribution and target-like contributions. In addition the 
experimental energy threshold is applied for detection of 
each isotope. The experimentally observed multiplicity 
distribution for most particles, including Z=l and 2, is 
well reproduced by the Rc= 5 calculation. When i?c=1.5 
is applied, the multiplicity is significantly overestimated 
for Z=l and 2 and underestimated for Z > 4. In the 
case of i?c=5, the calculated multiplicities start to devi- 
ate from those of the experiments at Z > 13, indicating 
that the projectile like contribution becomes significant 
for these fragments. 

Using the same prescription used for the experimental 
data, Oc/T, aa/T and ap/T are evaluated from the fil- 
tered yields of isobars. In FiglHl ln[R(I-|-2,I,A)] for the 
isobars with I=-l and 1 is plotted as a function of A for 
both primary and secondary fragments, together with 
the experimental results. No notable difference is ob- 
served between these two set of the calculated fragments. 
The calculated results are also well fitted by Eq.([n]) 
and the values for (/x„ — Hp)/T and Uc/T are extracted 
from the distributions of the primary and secondary iso- 
bars. The extracted values are (/x„ — ^p)/r=0.40 and 
ac/T=0.18 for the primary isobars and (/i„— y^p)/T=0.47 
and ac/T=0.17 for those of the secondary( Similar results 
are obtained from the unfiltered yields). The values fol- 
lowing secondary decay should be compared to the exper- 
imental values of (/x„ — ^p)/T— 0.71 and ac/T=0.35. The 
calcula ted values of — fJ-p)/T both from the primary 
and secondary fragments are somewhat lower than that 
of the experiments, as seen in the figure. The Coulomb 
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FIG. 9: Extracted values of symmetry energy coefficient from 
the experiments (solid circles) and calculations from the secondary 
fragments for Rc=5 (circles). Squares and triangles show those ob- 
tained for primary fragments for _Rc=1.5 and Rc=5, respectively. 
All errors evaluated are smaller than the size of the symbols. 



coefficient for the calculations is about a half of the exper- 
imental value, which suggests that the fragments in the 
calculations are more expanded and/or deformed than 
those of the experiments. 

In FigEl Ua/T values evaluated using Eq.dH) with the 
values of — Hp) /T and a^/T extracted from the model 
calculations are plotted as a function of A. The average 
values of the experimental results from Fig|4] (solid cir- 
cles) may be compared with those extracted following the 
de-excitation step in the Gemini Code(open circles). Al- 
though the experimental values show larger fluctuations, 
the calculated values and the data are in good agree- 
ment in magnitude and trend. In the figure the values 
extracted from the primary fragment yields are also plot- 
ted for both Rc = 1.5 (squares) and 5 (triangles). The 
values for Rc=5 show a slowly increasing distribution as a 
function of A with Oq/T ^ 4 to 5 over most of the range. 
Those for Rc—1.5 show a similar trend, but the values 
are smaller. The rather flat distribution of aa/T values 
of the primary fragments for A < 30 is consistent with 
a see nario of fragment emission from a common source 
with a given density and temperature. The heavier frag- 
ments with A > 30 may result from different mechanisms, 
e.g., projectile fragmentation in the more peripheral col- 
lisions. 

The pairing effect is experimentally observed in Fig|6l 
In the study of the complex fragments production in the 
238u ^ ].ga,ction at 1 A GeV, Ricciardi et al. suggested 
that the observed even-odd oscillation in the fragment 
yield originates from the last chance particle decay of the 
excited fragments during their cooling down evaporation 

m 



35| . In order to verify this hypothesis in 
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our experiment, we also did a similar study, using AMD- 
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FIG. 10: Calculated In[R(H-2,I,A)] values are plotted for 1=0 in 
the upper panel and for 1=2 in the lower panel. Triangles show 
results for fragments in the ground states, solid circles are for those 
in excited states just before the last particle decay. Dotted lines 
are to guide the eye. 



Gemini calculations. Two sets of simulated events were 
prepared. In one of the set, all fragments were in the 
ground state. In the other set, the last chance particle 
decay was blocked and all fragments were in an excited 
state just before the last chance particle decay. In the ac- 
tual calculations, the second set was generated from the 
first set by adding the charge and mass of the last chance 
emitted particle to that of the partner IMF. In FigflUl 
the results calculated for the different fragment sets are 
shown for ln[R(2,0,A)] and ln[R(4,2,A)] values in Eqs.® 
and ([TU|) . in which the pairing effect reveals, if it is suffi- 
ciently large. Calculated values of ln[R(2,0,A)] are shown 
in the upper and those of ln[R(4,2,A)] are in the lower 
figure. As one can see in the upper figure for the case 
of isobars with 1=0 and 2, the values calculated for the 
fragments in the ground state show a clear even-odd os- 
cillation pattern and the patterns are completely washed 
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out for fragments observed just before the last particle de- 
cay. This is consistent with the experimental observation 
shown in the upper part of FiglHl Even the disappear- 
ance of oscillations above A=30 is well reproduced. The 
same feature is also observed for isobars with 1=2 and 4, 
but the oscillation pattern is more distorted. In this case 
the results for the fragments in the ground state show the 
oscillatory pattern up to A = 28, but the pattern is not 
as clear as those in the experiments, especially for A < 
20, which are shown in the lower part of Fig El The os- 
cillatory pattern is totally smeared out for the fragments 
at an excited state just before the last chance particle 
emission, as in the 1=0 case. These results are consistent 
with those reported by Ricciardi et al. in their experi- 
ments and simulations, in which an enhanced oscillation 
pattern is observed for 1=0 isotopes and the pattern be- 
comes less clear for other isotope distributions [13, [l^ ■ 
Therefore our analysis strongly supports their hypothesis 
that the experimentally observed pairing effects in Figl6] 
indeed originate from the last chance particle decay dur- 
ing the statistical cooling down process of the excited 
fragments. 



V. SUMMARY 

For a large series of heavy ion reactions, coefficients 
of Coulomb energy, symmetry energy and pairing energy 
in the form of Oi/T as a function of A have been stud- 
ied from analyses of the yield ratios of isobars obtained 
in experiments and from model calculations. The AMD 
and Gemini codes were used for the calculations. For 
the symmetry energy term, the extracted values from 
the experiments are in good agreement with those calcu- 
lated for the final fragments in the ground state. They 
increase from 5 to ~ 16 as the masses of the fragments 
increase from 9 to 37. These values are generally much 
larger than those extracted from the primary fragments 
observed in the AMD calculations. Over the same mass 
interval the primary fragment values range from 4 to 5. 
This is consistent with a picture in which the primary 
fragments originate from a common emitting source. A 
smaller coalescence radius and earlier sampling time for 
the fragment formation in the AMD calculation result s in 
slightly smaller values of 2 to 3 for the primary fragments, 
but does not have a strong effect on the extracted values 
for the final fragments, though the isotope yield distribu- 
tion is not well reproduced compared to the case of Rc 
= 5. Although the technique employed is quite differ- 
ent, our model results are quite similar to those observed 
by Ono where the ratio C(Z) —Csym/T^ was extracted 
from the calculated isotope distributions for the reactions 
40ca-H40ca, 48ca+48ca, eoca+eoc^^ ^nd ^^Fe+'^'^Fe at 



yields from the three different reactions. The compar- 
isons between the experimentally extracted results and 
those of the calculations indicate that the experimental 
determination of symmetry energy coefficients, aa/T, are 
significantly affected b y the secondary decay processes 
of the primary fragments. This modification is a common 
feature of dynamic transport model approaches 36|, [37 1 . 
Thus extraction of the density dependence of the sym- 
metry energy from fragment observables must be done 
with caution and with appropriate attention to the role 
of the secondary decay. The importance of these effects 
will vary according to the observables employed for ex- 
traction of the desired information. 

The pairing effect is clearly observed in the experi- 
ments. Comparisons to the calculations strongly sup- 
port the hypothesis, which is proposed by Ricciardi et 
al. 34,13^1, that the observed effect originates from the 



last chance particle decay during the statistical cooling 
down process of the excited fragments. 

The Coulomb coefficient in the form of Qc/T is also 
evaluated both in the experiment and calculations. The 
experimentally extracted value is ac/T = 0.35, whereas 
the calculated values are Uc/T = 0.17 for the final frag- 
ments and Uc/T = 0.18 for the primary fragments. These 
difference suggest that the calculated fragments are more 
deformed and/or expanded than those observed in the 
experiments. 
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35 A MeV [36|. In that paper the ratios are evaluated 
from the quadratic shapes of the isotope distributions 
for a given Z after proper normalization of the isotope 
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